Origin of magnetic moments in defective Ti02 single crystals 
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In this paper we show that ferromagnetism can be induced in pure TiC>2 single crystals by oxygen 
ion irradiation. By combining x-ray diffraction, Raman-scattering, and electron spin resonance 
spectroscopy, a defect complex, i.e. Ti 3+ ions on the substitutional sites accompanied by oxygen 
vacancies, has been identified in irradiated TiC>2. This kind of defect complex results in a local 
(TiOe-a:) stretching Raman mode. We elucidate that Ti 3+ ions with one unpaired 3d electron 
provide the local magnetic moments. 



Recently, ferromagnetism has been observed in non- 
magnetically doped, but defective oxides, including 



TiO 



1.2.3.4 



This kind of observation challenges the 



conventional understanding of ferromagnetism, which is 
rather due to spin-split states or bands. Thus, one fun- 
damental question must be answered: where are the mo- 
ments located? Intensive theoretical work has been per- 
formed to understand the ferromagnetism in defective 
oxides£££ In these papers, the triplet states of p-likc 
electrons, located at cation or oxygen vacancies, yield 
the local moments, leading to a kind of ferromagnetism 
without the involvement of 3d electrons. Experimentally 
the ferromagnetism in undopcd TiC>2 has been found to 
relate with oxygen vacancies (Oy)^! 3 ., however, its mech- 
anism remains unclear. It is worth to note that Ti 3+ ions 
with one 3d electron are usually generated in slightly re- 
duced TiC>2. When O is removed, the excess electrons 
are unpaired^. They can occupy the nearby localized Ti 
3d orbit and therefore convert Ti 4+ ions to Ti 3+ ions. In 
a reduced rutile TiO^llO) surface, such a defect com- 
plex, Ti 3+ -Ov, has been well studied by first-principles 
calculation s 9 ' 10 and experimentally by resonant photo- 
electron diffractionii. Therefore, experimental work is 
needed to clarify whether the magnetic moments in de- 
fective Ti02 is due to unpaired 3d electrons localized on 
Ti 3+ ions. 

Ion irradiation is a non-equilibrium and reproducible 
method of inducing defects. Energetic ions displace 
atoms from their equilibrium lattice sites, thus creating 
mainly vacancies and intcrstitials. The amount of de- 
fects can be controlled by the ion fluence and energy. 
In this paper, we irradiated rutile Ti02 single crystals 
with 2-MeV O ions, resulting in a projected range of 
1.52 fim and a longitudinal straggling of 0.16 /*m as cal- 
culated by SRIM code (The Stopping and Range of Ions 
in Matter As a result of this irradiation, the forma- 
tion of Ti/O vacancies/interstitials is expected^. We 
selected high-energy oxygen ions as projectiles to avoid 
the introduction of foreign elements. Moreover, from a 
ballistic point of view, the creation of oxygen vacancies 
is more efficient, e.g., by a factor of 1.5 larger than the 
Ti-vacancy creation. From SRIM calculations it is also 
evident that, at the given energy, the maximum atomic 



concentration of the implanted oxygen ions is by a factor 
of 500 smaller than the concentration of oxygen recoils. 
For the region of maximum defect creation, i.e., at the 
end of the ion range, those projectiles play no chemical 
role. Thus, the desired local oxygen deficiency can be 
obtained by the proposed experimental procedure. The 
selection of TiC>2 single crystals instead of thin films is 
expected to exclude any effect of possible magnetic impu- 
rities in the underlying substrates^ 3 ., and of the interfaces 
between films and substrates^. 

Commercial rutile Ti02 single crystals were irradi- 
ated with 2 MeV O ions at room temperature with 
fluences from 1 x 10 15 to 5 x 10 16 cm -2 . All sam- 
ples are from the same purchase charge with an equal 
size of 5x5x0.5 mm 3 , corresponding to a mass of 53 
mg. The samples thereafter were named as 1E15, 
5E16. All samples were investigated using SQUID 
magnetometry (Quantum design MPMS-7XL), X-ray 
diffraction (XRD, Siemens D5005), Raman spectroscopy 
(LabramHR, Jobin-Yvon-Horiba) and electron spin res- 
onance (ESR, Bruker ELEXSYS E500 at 9.4 GHz). 
SQUID measurements indicate that the virgin Ti02 sin- 
gle crystals are weak paramagnets with a susceptibility 
of 7.7xl0" 6 cmu/T-g. 

Fig. QJa) shows the measured magnetic moment of 
Ti02 at 300 K before and after irradiation with a fluence 
of 5 x 10 15 cm -2 . After irradiation the sample shows 
a ferromagnetic hysteresis added to the linear param- 
agnetic background. The saturation magnetic moment 
is in the order of 10~ 5 emu, well above the sensitiv- 
ity limit of SQUID magnetometry. The inset shows the 
low-field part in an enlarged scale. The coercive field of 
around 10 mT can be clearly resolved. Fig. QJb) shows 
the irradiation-fluence dependent magnetic moment of 
Ti02 single crystals after subtracting the linear paramag- 
netic background. The magnetic moment first increases 
drastically with increasing fluence, but for the fluences 
above 5 x 10 15 cm~ 2 it decreases significantly again. For 
the 5E15 sample, with the largest magnetic moment, we 
also measured the hysteresis loops at different temper- 
atures (not shown). The loops are weakly temperature 
dependent. This is also a feature of the reported defect- 
induced fcrromagnetismi^. 
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FIG. 1: (Color online) (a) Magnetic moments measured at 
300 K as a function of magnetic field before and after irradi- 
ation. Inset: the low-field part of the loop for sample 5E15. 
(b) Magnetic moments at 300 K for samples with an equal 
sample-size and different ion fluences. The linear paramag- 
netic background has been subtracted. 



The structural evolution upon ion irradiation was as- 
sessed by X-ray diffraction (XRD) and Raman-scattering 
measurements. Fig. [2ja) shows the XRD patterns close 
to the TiO2(330) peak before and after irradiation at dif- 
ferent fluences. One remarkable characteristic is that a 
new peak appears at the left side of the (330) peak af- 
ter irradiation. This new peak is induced by a strained 
TiC>2 layer with a larger lattice spacing compared to that 
of the virgin TiC>2. With increasing fluence, this peak 
shifts towards smaller 2f? angles and arrives at a mini- 
mum at a fluence of 7.5 x 10 15 cm -2 . The largest strain is 
around 0.50%, however the peak becomes weaker. With 
further increasing fluence, the peak finally disappears and 
develops into a gradual shoulder, which means the loss 
of crystalline long-range ordering of the irradiated layer. 
This phenomenon is usually due to the fact that with 
increasing density point defects accumulate and develop 
into extended defects, e.g., dislocations. The lattice ex- 
pansion upon ion irradiation is a common feature in semi- 
conductors, e.g., in GaNi£. 

The effect of oxygen implantation on the structural 
properties of TiC>2 was further studied by Raman spec- 
troscopy. The experiments were carried out in a micro 
Raman 180° backscattcring geometry. The scattering in- 
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FIG. 2: (Color online) (a) XRD pattern for O-irradiated 
Ti02- The peaks (marked by arrows) in the left side of the 
TiO2(330) peak evidence the irradiation-generated strain in 
Ti02- (b) Raman spectra of Ti02 single crystals after oxygen 
ion irradiation with different fluences. The spectra are shifted 
vertically for better visualization. A clear mode at the range 
of Bi g (as indicated by the dashed arrow) has been observed 
after irradiation. 



tensities were measured using vertical orientation of the 
crystal [110]-axis relative to the electric field vector of the 
incoming laser. The spectrum of the virgin sample [Fig. 
HJb)] shows four distinct lines which are assigned to the 
B lg (143 cm" 1 ), E g (447 cm- 1 ), and A lg (612 cm- 1 ) fun- 
damental modes and to a higher order (HO) band around 
234 cm -1 of rutile TiO^- The fourth Raman-active fun- 
damental mode with Big symmetry appears as a weak 
shoulder at 820 cm -1 . Neither the fingerprint lines from 
the anatase Ti02 at 144 cm -1 nor those of Ti203 at 269 
cm- 1 and at 347 cm -1 have been observed in the irradi- 
ated sample o 18 ' 19 i 20 . Already at the smallest fluence, the 
irradiation induces i) the complete disappearance of B\ g 
and ii) significant changes in the range of the Bi g mode, 
where instead of the shoulder at 820 cm" 1 a clearly re- 
solved line appears at 834 cm -1 [indicated by the dashed 
arrows in Fig. HJb)]. The B\ g mode is very sensitive to 
the long-range order of Ti02 crystals. It softens when 
applying high pressure and looses intensity much faster 
than the other fundamental modes 2 ^,. Raman spectra of 
micro- and nanocrystalline Ti02 of 5 nm to 2 /im aver- 
age crystallite size did not show the B\ g mode as reported 
recently^. An important conclusion of the complete ab- 
sence of B\ g is that the information depth of the Raman 
spectra is limited to the irradiated part of the sample. 
The integrated intensity of the line at 834 cm- 1 is ~ 5.5 
times larger than that of the corresponding shoulder for 
a virgin Ti02- In virgin Ti02 the Bi g mode represents 
an asymmetric (Ti-O) stretching vibration of the (TiOg)- 
octahedra (Ti ions are surrounded by six O ions)2I. It is 
therefore sensitive to the short-range order and the bond 
strength in the crystal, i.e., the local environment of the 
Ti ions. Since a reduction of the crystallite size usually 
causes red-shifts of vibrational modes, the observed blue- 
shift by 14 cm -1 of the B^ g mode cannot be explained 
by phonon confinement. Furthermore, the line can be 
neither assigned to an IR-active nor to a silent mode of 
rutile Ti02 possibily activated due to a crystal-symmetry 
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FIG. 3: ESR signals observed in sample 5E15 at 5 K. The 
applied field is along the [100] axis. The main ESR peak 
accompanied by a hyperfine interaction pattern has a g-factor 
of 1.9408. 



reduction 2 ^. The line at 834 cm" 1 indicates a stronger 
bond between the Ti and the O atoms than in the virgin 
rutile. One possible reason is the reduced coordination 
number of Ti ions resulting in TiOe- x instead of TiOg, 
since a lower coordination number leads to bond contrac- 
tion and higher vibrational frequencies. This mechanism 
should in principle also affect the other Raman lines, par- 
ticularly the Ai g mode. However, a possible high-energy 
shoulder would be easily hidden under the large scatter- 
ing intensity of the main components. The localized ori- 
gin of the mode at 834 enr 1 is also supported by the de- 
pendence of the Raman line parameters on implantation 
fiuences. While the Ai g and E g lines are broadened with 
increasing oxygen fluence due to irradiation-induced de- 
fects, the local (TiOg-x) stretching mode shows a slight 
line narrowing. Therefore, we conclude that oxygen irra- 
diation results in the local distortion of the overall TiC>6 
coordination, which gives rise to a local TiOe- x stretch- 
ing mode in the range of the B^g mode. Such a distorted 
TiOe-x octahedral has also been considered in a recent 
experimental and molecular dynamics simulation study 
of ion irradiation damage in TiOz^ 3 -. 

In order to understand the origin of the local moments, 
ESR measurements were performed on the virgin and two 
selected irradiated samples: 5E15 (with the largest mag- 
netic moment), and 5E16 (with the largest irradiation 
fluence). ESR directly probes the interaction between 
unpaired electrons and the magnetic field. In a broad 
applied field range from 50 mT to 550 mT along the 
TiO2[100] axis and with temperatures down to 5 K, we 
only detected one ESR peak at a g-factor of 1.9408 ac- 
companied by a complicated hyperfine interaction pat- 
tern for the two irradiated samples. Fig. [3] shows the 
featured field-range of the measured ESR spectrum at 5 
K for sample 5E15. Note that electrons trapped by va- 
cancies have a <?-factor slightly larger than 2 and exhibit 
no hyperfine feature a 24 i 25 . The lineshape (see Fig. 2] for 
details) and the g-factor are typical features in the so- 
called A spectrum due to Ti 3+ ion o 26 ' 27 . Now we come 
to the question: how the Ti 3+ ions are generated after 
irradiation? Among the defects created by irradiation, 
only O vacancies and Ti interstitials can contribute ex- 
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FIG. 4: (Color online) (a) Angular-dependent resonance field 
of the ESR spectrum: it is split into two lines upon field 
orientation changing from the [100] to the [110] axis, (b)- 
(e) Comparison between simulated (solid lines) and measured 
(open circles) ESR data, (b) and (d): With field along the 
[110] direction without (b) and with (d) considering hyperfine 
interaction, (c) and (e): With field along the [100] direction 
without (c) and with (e) considering hyperfine interaction. 



cess electron o 8 ' 28 . There are three factors favoring the 
formation of Ti 3+ -Oy defect complexes: (1) As the near- 
est neighbors O vacancies are more close to substitutional 
Ti ions, (2) the number of O vacancies is 1.5 times larger 
than that of Ti interstitial*^ and (3) energetically the 
formation of O vacancies is slightly favored over the for- 
mation of Ti interstials 8 -. Therefore we can conclude that 
each of irradiation-induced oxygen vacancies leaves two 
electrons and converts the neig hboring Ti 4+ to Ti 3+ ions, 
resulting in Ti 3+ -Oy defect complexes. Such a convert- 
ing in the reduced TiO2(H0) surface has been confirmed 
recently by both theoretical and experimental work o 10 ' 11 . 
It is important to note that the virgin sample shows no 
detectable ESR line down to 5 K (not shown) since the 
virgin TiC>2 contains no unpaired electrons. 

The site occupation of the Ti 3+ ions, namely substi- 
tutional or interstitial, is determined by the angular de- 
pendence of the A spectrum. The crystalline symmetry 
results in a two-fold splitting of the Ti 3+ ESR lines at 
substitutional sites, and a four-fold splitting at intersti- 
tial sites, when the field orientation is changed from the 
[100] to the [110] axis^l. As shown in Fig. 0Ja), we can 
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conclude that the Ti 3+ ions in our samples occupy sub- 
stitutional sites. We performed a simulation of the ESR 
spectra with this assumption as shown in Fig. IHb)-(e). 
For fields along [110] and [100] and with the anisotropic g- 
factors [ga,, g y , g z }= [1.9408, 1.9738, 1.9766], we obtained 
the spectra as shown in Fig. B|b)-(d). Without consider- 
ing the hypcrfinc interaction, the simulated spectra agree 
reasonably well with the experimental data in g-factor, 
lincshape, and splitting with field rotation. When taking 
the hyperfine interaction into account, namely the nu- 
clear spin of 47 Ti (7 = 5/2) and 49 Ti (I = 7/2), the sim- 
ulations reproduce the experimental results even much 
better. 

The presence of Ti 3+ at substitutional sites also ex- 
plains the XRD results considering the fact that substi- 
tutional Ti 3+ ions are larger than Ti + ions resulting in 
the observed lattice expansion. However, the other point 
defects, i.e. Ti/O interstitials can also result in the lattice 
expansion. Combining the findings of Raman and ESR 
measurements, we conclude the following. Among the 
various defects created by oxygen ion irradiation, Ti 3+ 
ions on substitutional sites are accompanied by oxygen 
vacancies (Cv) forming defect complexes, which result 



in a local (TiOe- x ) stretching Raman mode. When the 
ion fluence is very large, the irradiated layer loses its 
long-range crystalline ordering, such that the fundamen- 
tal mode of Ti02 is broadened. 

Summarizing our results, we conclude that (1) a Ti 3+ - 
Cv defect complex is generated due to ion irradiation. 
The Ti 3+ ions provide local 3d moments which are de- 
cisively related with the observed ferromagnetism; (2) 
There is an optimum fluence to produce enough Ti + 
ions while keeping the crystalline ordering. Since ESR is 
not capable to determine the concentration of Ti 3+ ions, 
it is difficult to get a quantitative insight of the magnetic 
coupling mechanism. Despite such an unfolded question, 
our research suggests that the combination of ion irra- 
diation and ESR techniques can be a general approach 
to verify the predicted defect-induced ferromagnetism in 
oxides&2£. The simultaneous observation of ferromag- 
netism and ESR signal with a g-factor slightly above 2 
would be the criterion. 

The authors (S.Z. and H.S.) acknowledge financial 
funding from the Bundesministcrium fur Bildung und 
Forschung (FKZ03N8708). 



* Electronic address: S.Zhou@fzd.de 

1 J. M. D. Coey, M. Venkatesan, P. Stamenov, C. B. Fitzger- 
ald, and L. S. Dorneles, Phys. Rev. B 72, 024450 (2005). 

2 S. D. Yoon, Y. Chen, A. Yang, T. L. Goodrich, X. Zuo, 

D. A. Arena, K. Ziemer, C. Vittoria, and V. G. Harris, J. 
Phys.: Conden. Matter 18, L355 (2006). 

3 N. H. Hong, J. Sakai, N. Poirot, and V. Brize, Phys. Rev. 
B 73, 132404 (2006). 

4 S. Zhou, K. Potzger, G. Talut, H. Reuther, K. Kuepper, 
J. Grenzer, Q. Xu, A. Miicklich, M. Helm, J. Fassbender, 

E. Arenholz, J. Phys. D: Appl. Phys. 41, 105011 (2008). 

5 I. S. Elfimov, S. Yunoki, and G. A. Sawatzky, Phys. Rev. 
Lett. 89, 216403 (2002). 

6 J. Osorio-Guillen, S. Lany, S. V. Barabash, and A. Zunger, 
Phys. Rev. Lett. 96, 107203 (2006). 

7 T. Chanier, I. Opahle, M. Sargolzaei, R. Hayn, and 
M. Lannoo, Phys. Rev. Lett. 100, 026405 (2008). 

8 G. U. Oertzen and A. R. Gerson, J. Phys. Chem. Solids 
68, 324 (2007). 

9 P. J. D. Lindan, N. M. Harrison, M. J. Gillan, and J. A. 
White, Phys. Rev. B 55, 15919 (1997). 

I G. U. Oertzen and A. R. Gerson, Int. J. Quantum Chem 
106, 2054 (2006). 

II P. Kriiger, S. Bourgeois, B. Domenichini, H. Magnan, 
D. Chandesris, P. L. Fevre, A. M. Flank, J. Jupille, L. Flo- 
reano, A. Cossaro, A. Verdini, A. Morgante, Phys. Rev. 
Lett. 100, 055501 (2008). 

12 J. Ziegler, J. Biersack, and U. Littmark, The stopping and 
range of ions in matter (Pergamon, New York, 1985). 

13 F. Golmar, A. M. M. Navarro, C. E. R. Torres, F. H. 
Sanchez, F. D. Saccone, P. C. dos Santos Claro, G. A. 
Bem'tez, and P. L. Schilardi, Appl. Phys. Lett. 92, 262503 
(2008). 

14 A. Brinkman, M. Huijben, M. Van Zalk, J. Huijben, 



U. Zeitler, J. C. Maan, W. G. Van der Wiel, G. Rijnders, 
D. H. A. Blank, and H. Hilgenkamp, Nat. Mater. 6, 493 
(2007). 

15 G. Kopnov, Z. Vager, and R. Naaman, Adv. Mater. 19, 
925 (2007). 

16 C. Ronning, E. P. Carlsonb, and R. F. Davis, Phys. Rep. 
351, 349 (2001). 

17 S. P. S. Porto, P. A. Fleury, and T. C. Damen, Phys. Rev. 
154, 522 (1967). 

18 T. Ohsaka, F. Izumi, and Y. Fujiki, J. Raman Spectrosc. 
7, 321 (1978). 

19 A. Mooradian and P. M. Raccah, Phys. Rev. B 3, 4253 
(1971). 

20 G. A. Samara and P. S. Peercy, Phys. Rev. B 7, 1131 
(1973). 

21 V. Swamy, B. C. Muddle, and Q. Dai, Appl. Phys. Lett. 
89, 163118 (2006). 

22 V. Maroni, J. Phys. Chem. Sol. 49, 307 (1988). 

23 G. R. Lumpkin, K. L. Smith, M. G. Blackford, B. S. 
Thomas, K. R. Whittle, N. A. Marks, and N. J. Zaluzec, 
Phys. Rev. B 77, 214201 (2008). 

24 M. Sterrer, E. Fischbach, T. Risse, and H.-J. Freund, Phys. 
Rev. Lett. 94, 186101 (2005). 

25 S. V. Chong, K. Kadowaki, J. Xia, and H. Idriss, Appl. 
Phys. Lett. 92, 232502 (2008). 

26 P. F. Chester, J. Appl. Phys. 32, 2233 (1961). 

27 E. Yamaka and R. G. Barnes, Phys. Rev. 135, A144 
(1964). 

28 J. He, R. K. Behera, M. W. Finnis, X. Li, E. C. Dickey, 
S. R. Phillpot, and S. B. Sinnott, Acta Mater. 55, 4325 
(2007). 

29 P. Dev, Y. Xue, and P. Zhang, Phys. Rev. Lett. 100, 
117204 (2008). 



